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Abstract 

Collagen is the most abundant extracellular-network-forming protein in animal biology and is important in 
both natural and artificial tissues, where it serves as a material of great mechanical versatility. This versatility 
arises from its almost unique ability to remodel under applied loads into anisotropic and inhomogeneous 
structures. To explore the origins of this property, we develop a set of analysis tools and a novel experimental 
setup that probes the mechanical response of fibrous networks in a geometry that mimics a typical deformation 
profile imposed by cells in vivo. We observe strong fiber alignment and densification as a function of applied 
strain for both uncrosslinked and crosslinked collagenous networks. This alignment is found to be irreversibly 
imprinted in uncrosslinked collagen networks, suggesting a simple mechanism for tissue organization at the 
microscale. However, crosslinked networks display similar fiber alignment and the same geometrical properties 
as uncrosslinked gels, but with full reversibility. Plasticity is therefore not required to align fibers. On the 
contrary, our data show that this effect is part of the fundamental non-linear properties of fibrous biological 
networks. 
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Introduction 

Fiber networks, which arise in a range of natural and 
technological situations, are prime candidates for a 
wide spectrum of applications requiring tunable me- 
chanical, transport and chemical properties [1]. In 
nature, these networks self-assemble to generate im- 
portant structural and functional elements at vari- 
ous length scales: actin, intermediate filaments and 
microtubules are the main components of the cy- 
toskeleton [2]; spectrin confers versatile qualities to 
red blood cell membranes [3]; fibrin is an essential 
element in hemostasis [4]; collagen is the main com- 
ponent of the extracellular matrix (ECM) in the an- 
imal kingdom [5] and cellulose is used by plants to 
build cell walls [6]. 

The mechanical function of biological fiber net- 
works is essentially two-fold: (i) at the subcellular 
(actin, spectrin) and supracellular (collagen, fibrin) 
scales, the material offers little resistance and high 
sensitivity to small deformations, allowing it to be 
easily remodeled locally ; (ii) at larger strains it stiff- 
ens strongly to ensure cell and tissue integrity [7]. 



The non-linear stiffening, while observed in many bi- 
ological systems [8,9], is not fully understood yet, 
with theories focusing on one of two broad mecha- 
nisms: (i) microstructural nonlinearities of individ- 
ual filaments [8], and (ii) collective non-affine de- 
formations of multiple filaments [10, 11]. To un- 
ravel the relative importance of these mechanisms, 
a range of experimental tools have been developed 
to quantify the network's mechanical non-linearity 
in systematic ways and relate the material micro- 
structure (network density and morphology, fiber be- 
havior) to the mesoscopic stress-strain laws. These 
tools fall into two broad categories: simple shear in 
cone-plate or parallel plate geometries, and uniax- 
ial/biaxial stretch. 

Simple shear deformations are commonly used to 
study purified protein networks. This technique re- 
quires low sample volumes and provides a consis- 
tent set of experimental tools and generic protocols 
to probe the visco-elastic properties of soft gels in 
both the small-strain (linear) and large-strain (non- 
linear) regimes, and in addition, normal stresses can 
be measured. Recent data collected by Janmey et 
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al. [9] show in particular that sheared biopolymers 
exert negative normal forces, a fact that is in con- 
tradiction with the hyperelastic behavior of other 
well studied elastomers. The broad availability of 
experimental data in that geometry has encouraged 
a large number of related theoretical and numerical 
studies [12-14], focused primarily on the linear re- 
sponse of the material. However, since simple shear 
rheology assumes that the material undergoes purely 
isochoric deformations in the limit of small strains, 
it only allows for partial exploration of material be- 
havior. In particular, these experiments do not allow 
one to study completely the non-linear regime (strain 
typically larger than 10%) that is most relevant in 
many biological situations (single cell or tissue de- 
formation). And furthermore, it does not allow for a 
probe of the dilatational rheology of the networks. 

In contrast, at mesoscopic scales, uniaxial and bi- 
axial testing are most common for tissue mechan- 
ical characterization [7, 15-17] and have been used 
to study reconstituted collagen networks [18-22], the 
simplest tissue equivalents [23]. In contrast with sim- 
ple shear, uniaxial stretch generically leads to non- 
isochoric deformations, and hence allows one to mea- 
sure quantities such as the material Poisson ratio 
which can have values as large as v « 3 for strongly 
deformed collagen gels. These values arise in highly 
anisotropic materials, as reported for instance for 
solid foams [24], and it is somewhat surprising to 
see similar behavior in in vitro collagen gels which 
display little or no anisotropy in their undeformed 
state. 

To understand this, we recall that early studies [25] 
on cell/matrix interactions show that cells or groups 
of cells tend to generate tensile forces on the extracel- 
lular environment . When cell colonies were plated 
on fibrous materials such as collagenous gels, Har- 
ris and Stopak reported the formation of anisotropic 
and denser regions connecting these cellular assem- 
blies, and showed that the matrix structure has a 
strong influence on cell motility. Although these ob- 
servations are well accepted, little is known about 
the mechanical response of a fibrous matrix subject 
to an internal local strain. Neither of the mechanical 
characterizations described previously focus on how 
deformation changes the microstructure at the fiber 
scale, an issue of particular importance in the large 
strain regime, that is all too easy to observe (Figure 

In this paper, we use collagen type I gels as a 
model system to address this question and shed light 



on the morphological evolution of both the fiber [J 
and the network on an externally imposed stretch- 
ing strain. Collagen is a convenient biomaterial for 
biomechanical studies for a number of reasons: a) it 
is readily available in large amounts, which makes 
it suitable for milliliter-size gels; b) in vitro reconsti- 
tuted networks have fibers that are easily identifiable 
using confocal microscopy; c) many of its proper- 
ties have been extensively studied [18,26-30]; d) the 
large diameter of the fibers (« 0.5/xra for collagen 
fibrils [21]) and the stability of the network [31,32] 
make it easy to handle and image over a range of 
spatial and temporal scales; e) fibrillogenesis is con- 
veniently controlled in vitro by pH, temperature and 
concentration [26,28]. 
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Figure 1: Collagen gel morphological changes induced 
by presence of cells. (A) Single U87 glioblastoma cell in 
a collagen network 10 hours after gel polymerization. bar=50 
fim. (B) Several U87 cells on the surface of a collagen gel 
10 hours after gel polymerization. bar=200 fim. (C) Two 
cell colonies embedded in a collagen matrix 48 hours after gel 
polymerization. bar=200 /zm. Fibers (artificial red color) are 
imaged through confocal reflectance; cell nuclei (green) are 
labeled with a GFP-histone heterodimer. 

We first verified the presence of cell-induced align- 
ments and densification with our experimental sys- 
tem. As shown in figure [IjA., an isolated human 
glioblastoma cell (see Methods) in a collagen network 
induces stress variations and modifies the network 
texture in its vicinity [33,34]. Several isolated cells 
on the surface of a collagen gel produce fiber align- 

1 Although the expression "collagen fiber" traditionally 
refers to large bundles of collagen fibrils, we will use the words 
"fiber" and "fibril" interchangeably in this paper to refer to 
the 0.5 micron diameter bundles. 
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ment and network densification along lines connect- 
ing individual cells (figure J^B). Following Stopak and 
Harris, we also observed fiber alignment on macro- 
scopic length scales when we introduce large cell as- 
semblies in the same extra-cellular environment (fig- 
ure [Tp). Since active matrix remodeling is restricted 
to the vicinity of living cells, such an effect can only 
be accounted for by the mechanical properties of the 
network. 

With these observations in mind, we employ a 
specific experimental approach and develop a set of 
tools to quantitatively study the coupling between 
strain and the morphology of fibrous networks in a 
range of strain and strain rates that are typical of 
many biomechanical situations. Experiments on cell 
colonies suggest that such a process can be conve- 
niently studied at the millimeter scale, and over a 
time-scale of a few hours. However, instead of using 
cells to deform the extra-cellular matrix, we use an 
external imposed displacement to stretch collagenous 
samples and monitor the gel response. In particular, 
this experiment allows us to probe a range of dynami- 
cal regimes and independently tune the biochemistry 
(crosslinking) to study the coupling of tensile strain 
to network density and fiber orientation in a con- 
trolled setting and investigate the origin and general- 
ity of these mechanical processes. This also allows us 
to address the outstanding question of the mechani- 
cal reversibility of these patterns in an extracellular 
environment. 

Methods 

Network synthesis 

In vitro collagen networks are prepared according to 
a previously described cell culture-compatible pro- 
tocol [35], with a final collagen concentration rang- 
ing from 0.5 to 4.0mg/mL. Solutions consist of 10% 
10X minimum essential medium (Invitrogen, Carls- 
bad, CA), 10% fetal bovine serum (JRH Biosciences, 
Lenexa, KS), 1% penicillin-streptomycin (Invitro- 
gen), bovine collagen diluted to desired concentra- 
tion (from 3.1mg/mL or 6.4mg/mL batch, Inamed 
Biomaterials, Fremont, CA), a few /iL of 1M sodium 
hydroxide (NaOH, Sigma, St. Louis, MO) to bring 
pH to neutral, 50mM sodium bicarbonate (NaHC03, 
Sigma) buffer and deionized water. 800/iL of solution 
are pipetted onto glass-bottom Petri dishes (MatTek, 
Ashland, MA) Samples then polymerize for 30-60 
minutes in a cell incubator at 37° C, 5% CO2. Af- 
ter polymerization, samples are 20mm in diameter 



and l-2mm in height. 

In addition to untreated in vitro collagen gels, we 
also prepare polymerized samples, to which we add 
glutaraldehyde (GA) - a common cell and tissue fix- 
ative. The effect of this is an increase in the overall 
stiffness of the gel by at least an order of magni- 
tude (from a few tens of Pascals to over lkPa for a 
lmg/mL gel, as evaluted using a cone-plate rheome- 
ter), without noticeably changing the structure of the 
collagen network (fiber width and length, mesh size). 
2mL of 4% v/v GA (Sigma, St. Louis, MO) in water 
is pipetted onto the sample, which is then incubated 
once more for at least 2 hours. It is subsequently 
rinsed twice with deionized water. Before use, all 
samples are immersed into 2mL of deionized water, 
to allow the gel to swell, and to reduce friction. 

Cell experiments 

U87 human glioblastoma cells are cultured as de- 
scribed in [35]. After passaging the cells, tissue 
equivalents are generated by diluting the cells to 
approximately 10 5 /mL in an unpolymerized colla- 
gen solution at 1.5mg/mL. After polymerization in 
a 37° C temperature- and humidity-controlled incu- 
bator, the spacing between individual cells, as seen 
in figures |TJl and[I]3, is on the order of 100/im. Cell 
colonies (or spheroids), with an estimated 10 3 cells, 
are generated with the hanging droplet method [36] 
and subsequently seeded in 500/iL of collagen solu- 
tion at 1.5mg/mL shortly before it polymerizes. 

Bulk rheology 

We use an AR-G2 (TA Instruments, New Castle, 
DE) rheometer with a 4°, 40mm cone-plate geom- 
etry with a 109/im gap. 1.2mL of collagen solution 
is pipetted onto the 37° C preheated bottom plate of 
the rheometer and the cone is lowered onto the sam- 
ple. We use a solvent trap to prevent the sample 
from drying during the measurement. During poly- 
merization, the increase in G' and G" is probed by 
continuously oscillating the sample at a fixed 0.5% 
strain amplitude and at a frequency of 0.2Hz. The 
oscillatory strain sweep is performed at the same fre- 
quency and temperature, after the gel has polymer- 
ized for 2-3hrs. The strain amplitude is increased 
logarithmically until the sample breaks. 

During oscillatory strain sweeps, we simultane- 
ously record the maximum stress and strain of the 
sample for each oscillatory cycle. To characterize the 
onset of stiffening from the stress-strain data, we de- 
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fine the critical strain as the value at which the stress 
(7 exceeds the product G' j by more than 10%, where 
Gq is the elastic modulus in the small-strain linear 
regime. 



Mechanical setup 

We place a polymerized collagen sample onto a glass 
cover slip and perforate it with two rough-ended 
lmm-diameter glass cylinders (capillaries) (see fig- 
ures [2]A.,B), which are gently pushed all the way 
to the glass bottom to prevent the collagen from 
slipping beneath them. Each cylinder is attached 
to two secondary transverse elastic capillary rods, 
themselves attached to two linear transducers (New- 
port, Irvine, CA) controlled by the ESP300 controller 
(Newport). The transverse capillaries act as springs 
that allow to maintain contact with the bottom cover 
slip of the dish with constant pressure. The tips, ini- 
tially 1cm from each other, can then be moved apart 
at speeds ranging from 0.125 to 12.5 fim/s, corre- 
sponding to strain rates of 2.5-10 -5 to 2.5-10 -3 per 
second; this range includes measured rates of cell- 
induced contraction [37]. The movement of the tips 
results in the local stretching of the gel sitting be- 
tween them. For imaging purposes, the whole me- 
chanical setup (motors and tips) is clamped to the 
microscope sample holder plate. 



Confocal imaging 

We use a Zeiss LSM 510 Meta (Carl Zeiss Microimag- 
ing Inc., Thornwood, NY) equipped with a 488nm 
Argon laser line and several photomutliplier tubes. 
We set the Meta channel of the microscope (which 
allows for selection of specific wavelengths) to detect 
wavelengths between 474 and 494nm to allow for the 
fact that we work in reflectance mode [38] , which has 
the significant advantage of avoiding the use of fluo- 
rescent dyes in the samples. A 60X 1.2-NA Olympus 
water immersion objective (Olympus America Inc., 
Center Valley, PA) is mounted onto the microscope. 

While deforming the sample, we acquire timelapse 
2D confocal images at various heights between 50 and 
150/im from the bottom surface. Tracking in-plane 
deformations in multiple slices improves the statistics 
of our analysis; moreover, tracking out-of-plane mo- 
tion via image correlation allows us to estimate the 
vertical deformation of the sample (figure J2p). The 
timelapse interval is 10s, which corresponds to 0.25% 
imposed deformation at the typical strain rate. 
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Figure 2: Mechanical setup and sample imaging. (A) 

Side and (B) top views of the mechanical setup used to de- 
form the network; the collagen gel has a pancake-like shape, 
typically 1mm in thickness and 20mm in diameter. As defined 
in our experiments, the stretch axis is x. Drawn to scale, the 
two squares represent the flelds-of-view of the wide-field flu- 
orescence images (5x) and confocal reflectance images (60x). 
(C) Correlation of multiple slices over time gives an estimate 
of the interslice distances, and hence vertical strain. (D) Col- 
lagen network (blue) obtained with confocal reflectance. Flu- 
orescent tracers (pink) are embedded in the network. Scale 
bar 20/ira. (E) Wide-field fluorescence image of the embedded 
beads. Scale bar 500/ira. 



Image segmentation and fiber detection 

Fiber orientation is calculated for each 2D confocal 
slice by several image processing steps: the raw im- 
ages are filtered using a 2D Gaussian blur and subse- 
quently thresholded so that at 0% stretch, 10% of the 
pixels are above that threshold; this threshold value 
is applied to all subsequent images of the same ex- 
periment. A circular window of diameter 15px moves 
across the thresholded image, and the 2nd order mo- 
ment tensor M, defined below, is locally calculated 
using binary pixel weights: below the threshold level, 
the pixel weight is 0, above that level, A^ is 
1. Typically, the moment tensor quantifies the spa- 
tial distribution of weight around a center of mass 
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and its eigenvalues give an indication as to whether 
weights are distributed isotropically around the cen- 
ter or not. Similarly here, M quantifies the distribu- 
tion of pixel intensities around the intensity- weighted 
center of mass, with and the pixel coordinates 
in the local circular window: 

M \ IZAijiXij-x) 2 E^(x iJ -x)(y ii -y)| 

[EA t:j (X ZJ -X)(Y ZJ -Y) ^AijiYij-Y) 2 \ 1 ' 

Max {Eigenvalues(M)} 
; Min {Eigenvalues(M)} ^ ^ 

The ratio r characterizes the aspect ratio of bina- 
rized image fragment enclosed in the sliding window 
and can be used to detect fibers: a single fiber pass- 
ing through the middle of the window, with a diam- 
eter less than half the window size will yield a high 
aspect ratio; a single fiber in the corner of the win- 
dow or multiple fibers in the window will give a low 
aspect ratio. We only keep the regions with r > 2 
for which the eigenvector corresponding to the higher 
eigenvalue of M provides the local fiber orientation 

</>• 

The choice of window size (15 pixels) is a com- 
promise between increasing angular resolution at low 
strains and avoiding multiple fibers in the window. 
As shown later on, the fiber density strongly in- 
creases at high strain, resulting in less accuracy in 
the orientation analysis. In practice, this sets the 
limits of the method. See Methods SI, Figure SI 
and Figure S2 for validation of this image processing 
algorithm. 

Orientational order parameter 

To quantify the network anisotropy as a function of 
the applied deformation, we calculate an orientation 
tensor Q and an associated nematic order parameter 
\i from the distribution of the fiber orientation cj> at 
each time step, defined by: 

Q = f (cos 2 ((/))) (cos sin 0) \ ^ 
y (cos sin 0) (sin 2 </>) J 

fj, = Max {Eigenvalues(2 ft - Id)} (4) 

where Id is the identity matrix, and < • > denotes 
spatial averaging over the domain of interest. The 
order parameter ranges from \i = for a uniform an- 
gular (isotropic) distribution to \i — 1 for a perfectly 
aligned system. Although this parameter, based on a 
2D image analysis, only characterizes the order in the 



horizontal plane, it does provide a suitable signature 
of the microstructure evolution, and in particular its 
non-linear behavior. 

Deformation field at the mesoscale 

We also characterize, at the mesoscopic millimeter 
scale, the strain field e xx , e yy and e zz induced by 
the imposed displacement of the glass tips. Vertical 
strain is estimated by following the individual dis- 
placement of confocal z-stacks. For each slice of a 
stack taken at a time t, this is done by calculating 
the correlation with the slices obtained at a neighbor- 
ing time t-\-At (see figure [2p). The height of the slice 
which provides the largest correlation value indicates 
the new location of the material layer and this infor- 
mation is used to calculate the vertical component of 
the strain. 

The deformation in the xy plane is measured by 
a PIV (particle imaging velocimetry) method. \\im 
diameter rhodamine Fluospheres (Invitrogen, Carls- 
bad, CA), at a volume ratio of 1:100000, are used 
as tracers to measure the deformation of the gel at 
millimeter length-scales. Most of the carboxy-coated 
particles stick to the network, as seen in figure [2j3. 
Using a 5x lens on a Zeiss wide-field microscope and 
focusing on the middle of the sample, we image the 
local density of these particles, which displays het- 
erogeneities as seen in figure [2^1. An image cross- 
correlation technique is then used to track these het- 
erogeneities at a scale of 10 — 50/im as the network 
is progressively deformed. To identify the local dis- 
placement of a mesoscopic region of the gel located 
at (x, y) from time t to t + At, we extract a domain of 
48x48 pixels surrounding (x, y) at t and look for the 
best matching region - i.e. the one that maximizes 
the cross-correlation function - in the image obtained 
at t + At. The cross-correlation function is defined 
as: 

PAB = ^2 _ A avg)(Bij ~ B avg)/ (cr A B ) (5) 

where A and B are the pixel intensity values associ- 
ated with the two regions of interest, A avg and B avg 
are the local average pixel intensities, and a a and 
ob are the standard deviations of intensity values of 
those regions. 

This tracking through cross-correlation process is 
iterated over time to extract the full trajectory of 
the material point and corresponds to a Lagrangian 
description of the material. For a grid deformation 
that is fairly homogeneous in the field of view of the 
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microscope (millimeter scale) , we use a least-squares 
planar fit of the nodal displacement to quantify the 
material deformation at the mesoscale. The devia- 
tion from the fit provides a measure of the error on 
the deformation field. 

These measurements allow us to extract a number 
of strain and stress characteristics. In particular, in 
our geometry, the normal stresses <j yy and <j zz are 
negligible as the gel is not attached on the lateral 
sides, and this allows us to estimate the incremental 
Poisson ratio, defined as: 



-de y y/de x 



(6) 



which characterizes the coupling of incremental de- 
formations in orthogonal directions. In two dimen- 
sions and in the small deformation limit, the area of 
the grid A(e xx ) is related to the Poisson ratio via the 
following relationship: 

1 dA _ 1 
A(0)de xx ~ L Vxy [V 
Most materials respond to tension with a slight in- 
crease in their area (volume in three dimensions), 
which for an isotropic material translates into the 
condition: 



< 1 



(8) 

< 



The analogous condition in 3 dimensions is v xy 
0.5. At large deformations for an isotropic material, 
the criterion is slightly more complex, but the critical 
strain, beyond which the change of area (or volume 
in 3D) with respect to elongational strain becomes 
negative, remains of the same order in practice. 

Results 

Rheological characterization of the samples 

The range of collagen concentrations we work with 
(0.5-4.0 mg/mL) display mesh sizes from 1-5 /im 
(measured through analysis of confocal reflectance 
slices as in [35]) and span over two orders of mag- 
nitude in shear modulus (see figure [3]A.) , with val- 
ues in close agreement with previously reported data 
[35,39,40]. The linear shear modulus G' (measured 
at small deformation) has a strong dependence on the 
concentration c, with a behavior well approximated 
by G' ~ c 3 over the probed range. The onset of 
nonlinear strain-stiffening typically occurs at strains 
of the order of 5% (see figure [3^3 and [28]) and has 
only a weak dependence on collagen concentration 
(see figure [3K inset). 
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Figure 3: Bulk rheology measurements. (A) Linear elas- 
tic modulus as a function of collagen concentration. A power- 
law of 3 is shown for comparison. Inset: the onset of non-linear 
mechanical behavior, as defined in Methods. (B) Elastic mod- 
ulus G' during oscillatory strain sweeps and as a function of 
collagen concentration. 



Fiber alignment and non-linear Poisson effect 

Figures [4]A. and [4^ illustrate the evolution of the net- 
work microstucture over the domain as it is stretched. 
At low strains (< 5%), no particular alignment can 
be observed; however, above this threshold both fiber 
alignment and network density increase. This is 
quantified in the figures[4p and[4p where we show, as 
a function of the applied strain, the probability dis- 
tribution of local in-plane fiber orientation and the 
resulting order parameter. 

For a prescribed displacement of the capillary tips, 
we also characterized the gel microstructure as we 
move away from the axis connecting the two cap- 
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Figure 4: Typical results of a sample stretching experiment at micro- and meso-scale. The two montages of 5 images 
each show, for two different lmg/mL samples and at various strains, (A) wide-field fluorescence images of beads embedded within 
the network - the super-imposed grid is the result of the tracking of the deformation field (scale bar 500/im); (B) direct imaging of 
the fibers through CRLSM (scale bar 50/xm). In inset, each corresponding cf) histogram, with angle values going from to tt. For 
the same samples depicted above: (C) represents the evolution of the orientation statistics; the color at each point corresponds 
to the relative count of fibers oriented along a specific direction at a given strain e xx . (D) shows the order parameter /i resulting 
from the data in (C); the curve beyond 15% stretch is grayed out due to the lack of confidence of the order parameter when the 
high value of the density prevents a proper detection of the fibers (see Methods). (E) gives the deformations e yy and e zz as a 
function of the local strain e xx . In inset, the incremental Poisson ratio v xy as a function of the imposed deformation. 



illaries. Both the alignment and fiber density (ex- 
pressed as fraction of pixels above a given thresh- 
old) decrease (figure [5]). This picture is, as expected, 
in direct agreement with the alignment pattern in- 
duced by cell colonies pulling on extra-cellular matrix 
shown in figure [I] where fiber alignment and density 
are maximum along the axis joining the colonies, and 
decay away from it. 

In figure [4^1 we show the average induced strain 
components e yy and e zz as a function of the exter- 
nally applied strain e xx . When e xx < 2% there is 
little contraction in the transverse direction so that 
the Poisson ratio v xy ~ 0. As the applied strain in- 
creases, the material first thins by contracting in the 
z direction when 2% < e xx < 5%, and only when 
£ X x > 5% does it also contract in the transverse y 
direction, with observed values of v xy as high as 5 
(figure ptl inset). The lag in response between these 



two directions can be attributed to the sample ge- 
ometry as well as a slight initial anisotropy of fiber 
orientations in the yz plane [21]; here we consider 
only the properties of the planar projection of the 
network. The large in-plane incremental Poisson ra- 
tio v xy quantifies the change in local fiber density and 
is consistent with the confocal observations of den- 
sification. In order to quantify and compare these 
changes with an independent measure of the geom- 
etry of deformation, we define a critical strain a C rit 
for which dA(e cr n) j de xx = 0, beyond which the areal 
strain (yellow frame on figure [4]A.) starts to decrease 
with the applied strain. We find that e cr i t ~ 5%, 
consistent with the critical strain observed for fiber 
alignment. 

The critical deformation, as defined above from the 
kinematic behavior in local stretching tests, can be 
directly compared with the strain associated with the 
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-2000 -1000 1000 2000 
distance from pulling axis [microns] 

Figure 5: Characterization of the spatial variations of 
collagen fiber alignment and densification. (A-C) Con- 
focal reflectance images of lmg/mL collagen network stretched 
up to 15%; images are located at (A) 0, (B) 0.5 and (C) 2 
millimeters from the stretching axis, bars = 50/ira. (D) Or- 
der parameter and density as a function of distance from the 
stretching axis, for the same sample. Letters A, B and C on 
figures (D) and (E) refer to corresponding images above. (E) 
Drawn to scale, locations of images (A), (B) and (C) with 
respect to stretching axis. 

mechanical stiffening measured in rheological experi- 
ments (see figure [3]). These two quantities, measured 
independently, show good correlations in their val- 
ues and trends. For unfixed collagen samples, crit- 
ical strain values, measured either from rheological 
meansurements or from the kinematics, range from a 
few percents at high collagen concentration to 15 % 
at low concentration (see figure|6]). The critical strain 
is therefore very weakly dependant on the concentra- 
tion, in particular compared with the variation of the 
elastic modulus at small deformation that varies over 
more than two orders of magnitude in the same con- 
centration range. G A- fixed samples, whose stiffness 
is estimated at an order of magnitude higher than 
their non- fixed counterparts, also show similar val- 
ues for e cr it. Taken together, these results show that 
the strains above which the gel behavior becomes 
non-linear as evidenced i) from the elastic modulus 
for a simple shear geometry and ii) from the Poisson 
ratio in the local stretching experiments are related 
with each other and only weakly sensitive to physical 
scales such as the actual value of the elastic modulus. 

Orientational ordering is an elastic effect 

Strain-induced alignment arises a priori from a com- 
bination of reversible elastic effects and irreversible 
inelastic effects. To disentangle these two contribu- 
tions, we apply repeated strain cycles to the sample. 
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Figure 6: Comparison of critical strains measured from 
bulk rheology and 2-point stretching. Comparison of 
critical strain measured from the PIV method (e C rit such that 
dA(e cr it)/de xx = 0, see Results) on locally stretched samples 
with the onset of strain stiffening (~y C rit such that a > l.lG f j, 
see Results) obtained from rheological measurements. 

All pure type-I collagen samples display very little re- 
versibility regardless of their concentration once the 
imposed strains exceed about 10%; the gel never re- 
covers its initial configuration once the capillary tips 
return to their initial location. In figures [7]A.,B we 
show the evolution of the fiber orientation histograms 
as the material is cyclically stretched up to strains of 
15%: fiber aligment is permanently imprinted (J7p- 
E). 

With the addition, pos£-polymerization, of GA to 
our samples (see Methods), we change collagen's ma- 
terial properties (elastic modulus, plasticity thresh- 
old) without changing the microstructure of the net- 
work. For small deformations, the initial response is 
similar to that of the unfixed sample, indicating that 
fiber alignment and the anomalous Poisson effect are 
only weakly sensitive to the fiber material proper- 
ties. Furthermore, during cycles of applied defor- 
mation - i.e. ramping up the applied strain to 15% 
and returning back to 0% at the same rate - GA- 
crosslinked samples exhibit near-reversibility at the 
microscale, with fiber images of successive cycles be- 
ing almost identical (figures [7}H-J). Figures^ and 
[7p quantify this reversibility in terms of the orien- 
tation histogram and the order parameter over mul- 
tiple cycles. Consistent with this behavior, we find 
(see figure [TjK) that the local deformation field at 
the mesoscale evolves along a reversible path for the 
same strain rate (2.5- 10 _4 /s). Taken together, these 
results show that fiber plasticity, though observed 
for pure collagen samples, is unimportant in deter- 
mining alignment at microscopic scales and the large 
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Figure 7: Characterization of the reversibility of alignment and densification in crosslinked and uncrosslinked 
collagen samples. Response to cycles of deformation for untreated collagen gels (A-E) and gels treated with glutaraldehyde 
(F-K). (A) Histograms of the fiber orientation as a function of the imposed strain for untreated sample; color at each point 
corresponds to the relative count of fibers oriented along a specific direction at a given strain e xx . (B) Resulting order parameter 
for untreated sample. (C-E) Confocal reflectance images of a lmg/mL collagen sample cycled 4 times up to 15% stretch and back 
to 0%, showing the extent of the reversibility at the microstructural level: (C) beginning of 1st stretch cycle; (D) middle of 1st 
stretch cycle; (E) end of 1st stretch cycle; bars = 50/ira. (F) Histograms of the fiber orientation as a function of imposed strain 
for glutaraldehyde-treated sample; color-coding as in (A). (G) Resulting order parameter and density for glutaraldehyde-treated 
sample. (H-J) Confocal images of a lmg/mL collagen samples with glutaraldehyde, cycled 4 times up to 15% stretch and back to 
0%: (H) beginning of 3rd stretch cycle; (I) middle of 3rd stretch cycle; (J) end of 3rd stretch cycle; bars = 25/ira. (K) Mesoscopic 
response to 15% strain cycles for various strain rates at a constant imposed strain amplitude. The bold line corresponds to the 
strain rate used for all other experiments. 



Poisson ratio at larger scales. That is, strain-induced 
alignment is primarily an elastic effect. 

The reversible behavior characterized above might 
then serve as a baseline to study more subtle effects 
(e.g. time-dependent and/or irreversible processes) 
that can be observed after many cycles or different 
strain rates. We see, in particular, a slight decrease 
in the amplitude of the alignment with the number 
of cycles, indicating that fiber plasticity still occurs, 
although on much larger timescales. In figure [7]K 
we show the effects of strain rate on the response of 
the system for a fixed strain amplitude. At lower 
strain rates the system does not recover completely 
after a full cycle. This offset in the response can 
be attributed to a slow creeping process occurring 
over a time-scale of a few hours. At larger strain 
rates, one expects to see dynamic effects related to 
viscous dissipation. Previous studies on the poroe- 
lasticity of collagen networks [18,20,41-43], have re- 
ported equilibration times ranging from a few sec- 



onds to a few hours, probably due to the diversity 
of geometries and setups used for all these measure- 
ments. We find that dynamic effects associated with 
higher strain rates induce an asymmetry in the re- 
sponse between loading and unloading, influencing 
the unloading curve more than the loading curve. 
This shows that the material responds with differ- 
ent time-scales in extension and compression; this, 
in turn, suggests that different physical processes are 
involved during the loading which is dominated by 
fiber stretching and unloading which is dominated 
by fiber bending. 

Discussion 

In this study, we used an experimental and com- 
putational approach to quantify the emergence of 
fiber alignment as a collagen sample is stretched. 
This behavior is consistent with observations of cell- 
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induced morphological changes in tissue equivalents 
and sheds light on a biologically relevant material 
non-linearity that arises from stress heterogeneities 
in fiber networks. Fiber alignment at the microscale 
results in tissue densification when boundary condi- 
tions allow it; this leads to high values of the Poisson 
ratio at large deformations, and is observed for all 
studied collagen concentrations, with or without ad- 
dition of glutaraldehyde. 

The reversibility of fiber alignment and gel densi- 
fication, seen in crosslinked collagen samples, show 
that these effects are primarily elastic. Experiments 
on a piece of synthetic felt [44] have demonstrated 
that geometry alone can account for such a behav- 
ior based on the generic non-linearity of individual 
fibers, stiff in extension, but soft in compression 
(bending/buckling). Under uniaxial loading, a ten- 
sile stress is necessarily balanced in the microstruc- 
ture by a compressive load on fibers normal to the 
stretch direction, leading to a collapse of the material 
in the normal direction and a strong enghancement 
of fiber alignment along the load direction, as ob- 
served here. This argument also explains the correla- 
tion between the moment where stress builds up (on- 
set of non-linearity in rheological measurements) and 
the critical strain associated with the Poisson effect. 
This is also consistent with the bulk rheology exper- 
iments performed by Janmey et al. [9] who report 
large negative normal stresses in all tested fibrous 
materials. The normal stress in volume-constrained 
geometries (as in simple shear flows) is the counter- 
part of high Poisson ratios in unconstrained tests 
such as the local stretch performed in the present 
study. 

Our findings seem to contrast with the experi- 
mental results of Tower et al. [19] on stretched col- 
lagen samples, where either irreversible alignment 
(pure collagen sample) or early fracture with negli- 
gible alignment (GA crosslinked) occurs, suggesting 
that fiber plasticity is a key player in the alignment 
process. The difference with the present study can 
be readily explained from a geometrical standpoint: 
Tower et al. used a geometry in which the width 
of the stretched region is comparable to its length; 
this might significantly constrain transverse motion 
of the material and prevent local volume changes to 
their full extent. By contrast, the 2-point stretch- 
ing device we use ensures that the material is free to 
move along the direction normal to the load. 

Our experiments have demonstrated the impor- 
tance of sample geometry and boundary conditions 
on the microstructure and mechanical response of 



reconstituted biopolymer gels. For functional tis- 
sues, it is known that mechanical properties are of- 
ten finely controlled by the texture of the underlying 
extracellular matrix as well [45]. Understanding the 
mechanisms leading to such organization is an impor- 
tant step in learning how it happens in the formation 
of natural tissues and for developing strategies to en- 
gineer suitable tissue equivalents. We have shown 
here that ECM texturization can be brought about 
simply by applying a deformation with a purely me- 
chanical device, without any intervention due to ac- 
tive modeling by cells. This external perturbation is 
applied pos£-polymerization, in contrast with previ- 
ous reports of fiber alignment induced by a flow [46] 
or a magnetic field [47, 48] during collagen polymer- 
ization. This provides direct support to the in vivo 
studies of pos£-polymerization collagen texturization 
in developing tissues [49]. However, the microscopic 
origin of the permanent texturization occurring in 
our uncrosslinked samples remains unclear. Fibril- 
fibril junctions are likely to be where plastic defor- 
mation occurs, allowing fibrils to slide with respect 
to one another and thus inducing irreversible changes 
in the network topology. One of the effects of glu- 
taraldehyde crosslinking is to strengthen these junc- 
tions and reduce the amount of plastic reorganization 
allowed at the network level. 

Whatever their origin, geometrical changes in the 
network structure are known to influence, in turn, 
cell behavior. Our mechanical setup allows for the 
dynamical control of network texture in a passive 
system, but clearly can be used to study how ac- 
tive cells respond to externally-induced anisotropy. 
Such experiments will provide more insights into spe- 
cific mechanisms of mechanotransduction and cell be- 
havior, which are crucial to processes such as mor- 
phogenesis, stem cell differentiation, metastasis and 
wound healing. More generally, the external con- 
trol of fiber alignment post polymerization not only 
offers a convenient way to design anisotropic tissue 
equivalents with collagen, but can also, when applied 
to other biopolymer systems, shed light on a range 
of analogous phenomena, such as actin gel contrac- 
tion [50] or platelet-fibrin interactions [51], where mi- 
croscopic agents interact through the network and 
lead to large scale evolution and reorganization of 
matter. 
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